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Abstract. Reaction-diffusion analytical modeling of predator-prey systems has shown that specialist natural ene- 
y^mies can slow, stop and even reverse pest invasions, assuming that the prey population displays a strong Alice 
effect in its growth. Few additional analytical results have been obtained for other spatially distributed predator- 
I prey systems, as traveling waves of non-monotonous systems are notoriously difficult to obtain. Traveling waves 
*^have indeed recently been shown to exist in predator-prey systems, but the direction of the wave, an essential item 
■^-nof information in the context of the control of biological invasions, is generally unknown. Preliminary numerical 
(“^ explorations have hinted that control by generalist predators might be possible for prey populations displaying 
■^logistic growth. We aimed to formalize the conditions in which spatial biological control can be achieved by gene- 
gralists, through an analytical approach based on reaction-diffusion equations. 

The population of the focal prey — the invader — is assumed to grow according to a logistic function. The 
OQpredator has a type II functional response and is present everywhere in the domain, at its carrying capacity, 
^on alternative hosts. Control, defined as the invader becoming extinct in the domain, may result from spatially 
^^independent demographic dynamics or from a spatial extinction wave. Using comparison principles, we obtain 
,__|Sufficient conditions for control and for invasion, based on scalar bistable partial differential equations (PDFs). 
^NjThe searching efficiency and functional response plateau of the predator are identified as the main parameters 
^^defining the parameter space for prey extinction and invasion. Numerical explorations are carried out in the region 
those control parameters space between the super- and subsolutions, in which no conclusion about controllabi- 
f^htv can be drawn on the basis of analytical solutions. 

in 

'. ^ The ability of generalist predators to control prey populations with logistic growth lies in the bistable dynamics 

the coupled system, rather than in the bistability of prey-only dynamics as observed for specialist predators 
^attacking prey populations displaying Allee effects. Analysis of the ordinary differential equations (DDEs) system 
^identifies parameter regions with monostable (extinction) and bistable (extinction or invasion) dynamics. By 


a 


contrast, analysis of the associated PDE system distinguishes different and additional regions of invasion and 
extinction. Depending on the relative positions of these different zones, four patterns of spatial dynamics can be 
identified : traveling waves of extinction and invasion, pulse waves of extinction and heterogeneous stationary 
positive solutions of the Turing type. As a consequence, prey control is predicted to be possible when space is 
considered in additional situations other than those identified without considering space. The reverse situation is 
also possible. None of these considerations apply to spatial predator-prey systems with specialist natural enemies. 
The consideration of space in predator-prey systems involving generalist predators with a parabolic functional 
response is thus crucial. 
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1 Introduction 


1.1 Modeling the biological control of invasive pests 


Biological invasions are a major contemporary problem 

(Pimentel |2011 

Gamier et al. ||2012 

Mistro 

et al. 

2012 

Potapov &: Rajakaryne 

2013 

Wang et al. 

2013 , [Savage Sz Renton 

2013 

) for which few 


solutions are available, all of which are very costly. The use of natural enemies for the biological control 


of invading insects is one of the most promising possibilities (Moffat et al. 2013, Li et al. 2014 Ye et 


al. 2014, Basnet &: Mukhopadhyay 2014). As invasion is essentially a spatial process, the potential of 


natural enemies to stop or even reverse an invasion is of particular interest. The fundamental analytical 
work of |Owen fc Lewisl2001 showed that specialist predators could potentially slow, stop or reverse the 
spread of invasive pests. The reversal of pest spread by specialist predation requires a strong Allee effect 
for the pest-only dynamics, defined as a negative growth rate for the prey population at low density. In 
the presence of a weak Allee effect, the predator can stop, but not reverse the wave of invasion. These 


conclusions have been confirmed in several other theoretical studies (Cai et al. 2014; Boukal et al. 


2007, Morozov Sz Petrovskii 2009). 


Generalist predators can also control prey effectively (Erbach et al. ]2014t [Chakraborty [2015] ). Their 
use could be promoted through conservation biological control programs without the need for exogenous 
specialist natural enemies. Unfortunately, the role of generalist predators in the spatial control of their 
prey has been much less studied than that of specialist predators, due to the intrinsic difficulties of having 
to work with a system of equations rather than with a single scalar equation. However, two important 
studies have been carried out in this area : the analytical and comprehensive study of Du & Shi 2007 


and the preliminary simulation study of |Fagan et al. |2002 Both used the same model structure as we do 
here, with logistic growth for both prey and predator populations, and a type II functional response for 


predators. The convergence of these models was strengthened further by the in-depth analysis of Magal 


|et al. |2008 in which space was not considered. It is difficult to use these models in a spatial context : the 
work of Du &: Shr||2007 cannot deal with invasion and traveling waves, because it deals with a bounded 
space. The numerical simulations of Fagan et al. 2002| are restricted to a few parameter values. They 
are, however, valuable, because they suggest conditions in which a generalist predator might be able to 
stop, and even reverse the invasion wave of a pest population displaying logistic growth. Fagan et al. also 
reported the results of field studies indicating that predators with diffusion coefficients higher than those 
of their prey are poor control organisms. The authors provided an explanation for this finding founded 
on logical arguments, but without a firm mathematical foundation. This result has been confirmed 
by a few numerical simulations including space, as reported by Magal et al. 2008 revealing a strong 


dependence of system dynamics on the relative rates of diffusion of the prey and the predator. It is 
thus important to take space into account, by considering diffusion coefficients of both predator and 


prey. This conclusion accounts for the interest of scientists in questions of this type (Lewis et al. 2013 


Hastings 2000; De Roos et al. 1991, 1998). 


There are therefore hopes that it might be possible to extend the conditions for the control of inva¬ 
sive prey organisms to (i) generalist predators and (ii) prey populations displaying growth patterns not 
dependent on the restrictive assumption of Allee effects. Such control approaches would have a major 
impact in the field, given the high degree of generalism obtained. The aim of this study was, therefore 
to formalize the conditions in which spatial biological control can be achieved by generalists, through 
an analytical approach based on traveling waves solutions of reaction-diffusion equations. 

Traveling wave solution describes a constant profile U moving through space at a speed c. Such waves 
are often observed in nonlinear reaction-diffusion systems modeling various phenomena. They are par¬ 
ticularly suitable for describing the propagation of invasive fronts. In systems modeling a single species, 
described by a scalar equation, this type of solution is very well understood (Fischer |l937 Kolmogorov 


et al. 1937 and Volpert et al. 1994 for a complete theory). Two particular classes of equations can 


be distinguished : monostable equations (like the Fisher-KPP equations) and bistable equations (often 
modeling the Allee effect). In monostable equations, there is a minimal wave speed c* such that, for 
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any c > c*, a wave solution with speed c exists. In bistable equations traveling waves exist for a unique 
speed c*. The sign of this speed c* distinguishes between invasion or extinction of prey, which is a key 
property for our purposes. 

For interactions of several species (described by a multidimensional system), the situation is much 
more complex. However, for some type of interaction, cooperation for instance, the system possesses 
a strong structural property, namely monotonicity. Essentially, this monotony makes it possible to use 
the comparison principle, which is always possible for one-dimensional systems, and the theory is then 
complete (see Volpert et al. 1994). Unfortunately, our system, and prey-predator systems in general, 


do not have such a monotonous structure. This method is then unsuitable for monotonous systems and 
only a few results have been published. One of the key reasons for this is as follows : when we search for 
traveling wave solutions for a system with N equations, we obtain a system of N second order ordinary 
equations that can be reformulated as a system of 2N first order ordinary equations. In the scalar case 
(N = 1) , it is therefore possible to study trajectories in a plane, available using classical tools for 
two-dimensional dynamical systems. For several species {N > 1), it is necessary to study trajectories in 
a 2A^-dimensional space, which may be very difficult. 

Hence, the first rigorous results demonstrating the existence of traveling waves in prey-predator 
systems were based on a generalization, to the fourth dimension, of the classical shooting method in the 

1984a|b). This approach has since been generalised (Huang et al. 


phase plane (Dunbar 


2003; Xu & 


Weng 2012). However, all these studies simply investigate the mere existence of traveling waves. They 
do not determine the direction of the wave or the global dynamics for general initial conditions. Other 


methods have recently been developed in similar models (Huang &: Weng 2013; Ducrot & Langlais 


2012), but they are subject to the same limitations. A last approach is to use the degree theory (see 


Giovangigli ||1990 Volpert et al. 1994) to obtain the existence of traveling waves. These homotopy 


e.g. 

methods may occasionally give some information on the speed c. Unfortunately, this needs additional 
estimates which are very difficult to obtain here. We therefore required another method. 


The analysis provided in Magal et al. 2008 gives conditions for preys’ control by predators, but 


this analysis was carried out largely without reference to space. Thus, we have extended the system of 
[Magal et al. |[2008 by adding spatial diffusion. We find that the conditions for control are very different 
from those for the system in which space is not considered. The conditions for prey extinction and 
invasion are discussed in terms of two essential parameters : the encouter rate E and the handling time 
h. Increasing E clearly increases predator pressure. Conversely, increasing h decreases predator pressure. 


The paper is organized as follows. In section we present the mathematical model and the main 
result of this work : theorem 2.1 describes invasion conditions for the ODE system and the theorems 


2.4 and 2.6 the invasion conditions for the PDE system. The mathematical results are completed by 


numerical simulations in section |3l The results are discussed in section [H The final section [5] is devoted 
to the mathematical proofs. 


2 Model and main results 

2.1 Mathematical model 

We analyze a system of partial differential equations for a prey population with logistic growth, and 
a generalist predator population with logistic growth on alternative prey in the absence of the invading 
host. The functional response is of Holling type H. The prey-predator interactions are modeled by the 
following partial differential equation system : 

dtu = DuAxU + riu 

< ^^v = D^A^v+ r2V (l - x G M, t G M+ (1) 

tt(x, 0) = Mo(a^)5 u(x,0) = l 
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with : 

u{t, x) = prey density at time t and at point x. 
Du = diffusion rate of prey 
ri = growth rate of prey 
Ki = carrying capacity of prey 
E = encounter rate 
7 = conversion efficiency 


v{t,x) = predator density at time t and at point x. 

Du = diffusion rate of predators 
r 2 = growth rate of predators 

K 2 = carrying capacity of predators in absence of focal prey 
h = handling time 

uq >0 the initial concentration of prey, 


Du, Du, ri,r 2 , Ki,K 2 , E, h and 7 are positive constant parameters. 


We carried out the following adimensionalization ; 

t' = Tit ; x' = ; u'{x', t') = u{t, x)jKi ; v'{x', t') = v{t, x)lK2 

d' = Du/Du ; r' = ra/n ; E' = EK2/ri ; h' = rihKi/K2 ; 7' = 1K1/K2 ; a 
Removing the sign ’ to simplify the notation, the system reads 


J dtu = l:^xU + u{l-u) - xGM,tGM+ 

I dtv = dA^v + r(v{l-v) + 

with the initial conditionsQ 


n(0,x) = Uo{x) G [0,1] ; 
u(0, x) = 1. 


lim uo{x) = 1 

03^ — 00 


lim ^ 0 ( 3 ;) = 0 

ir^+oo 


( 2 ) 

(3) 


2.2 Main results 

We distinguish two ways in which a predator can control the prey, one taking space into account 
and the other not considering this factor (mathematical definitions are provided in definition 2.2.1). 

- The spatially uniform extinction results exclusively from local demographic processes and is inde¬ 
pendent of space. 

- The extinction wave is due to both demographic and diffusive processes and may take various 
forms, from a traveling front to a pulse. 

Conversly, invasion is defined as prey survival and we distinguish two ways in which the prey can invade. 

- The spatially uniform invasion, which is independent of space. 

- The non-uniform invasion, described by various spatial dynamics, from Turing phenomena to 
invasion waves. 

Definition 2.2.1. Let { uq { x ), vq { x )) be an initial condition verifying (§ and {u{t,x),v{t,x)) be the 
corresponding solution of Q. 

- Extinction of prey oecurs if 

VxGM, lim u(t,x) = 0. 

t^-eoo 

- Prey extinction is uniform if it is uniform with respect to x G M, that is, if there exists a map 4>{t) 
verifying 

Vx G M, Vt > to) 0 < u{t, x) < (/{t) and lim (p{t) = 0. 

t^-l-oo 

- Prey extinction is non uniform if there is extinction but no uniform extinction. 

- Invasion of prey occurs if there is no extinction, that is if 


3x G M, limsupu(t, x) > 0 

t —^“|~00 


1. All our results remain true for various different initial conditions. The essential condition is that the solutions of the 
scalar systems we consider converge to traveling wave solutions. In particular, compact support may be allowed for uq. See 
Fife ||1979 for a detailed discussion. 
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2.2.1 Analysis of the associated ODE system 

If space is not taken into account, system Q may be rewritten as follows. 


|u = n(l-u)- 


Euv 

1+Ehu'^ 

Euv 


t e 


iv = r [v {1 - v) + ) . (4) 

< 1 < 

2008[). Indeed, it is clear that there are always three 


' l-\-Ehu J 

0 < n(0) = tto < 1 < u(0) = vq 

System @ is well understood ( |Magal et al 
trivial stationary states : (0,0) and (1,0), which are unstable and (0,1), which is asymptotically stable 
if, and only if, E > 1. Moreover, there are no more than three non-trivial positive steady states. We 
are interested principally in the case E > 1. In this case, there are either no or two stationary positive 
steady states. If the two steady states exist, denoted {u,v) and {u*,v*) with u < u* and v < v*, then 
{u,v) is always unstable and {u*,v*) is most often stable. In this case, there are two stable nonnegative 
solutions, (0,1) and {u*,v*) and the system is bistable. 

We are interested principally in the conditions for prey extinction. If E < 1, then (0,1) is unstable 
and no extinction occurs. We are therefore interested only in the case E > 1. Now, if E > 1, there are 
two possibilities. In the non bistable case, (0,1) is globally stable and there is extinction. In the bistable 
case, provided that u is initially small enough, say uq < fii for some 0 < /ri < 1, then u{t) —>■ 0 as 
t —>■ + 00 . Conversely, if u is initially large enough, say 0 < /i 2 < uq, then u{t) —>■ u* as t —>■ +oo. Thus, 
in this case, the outcome — extinction or invasion — depends on the initial conditions. The following 
result provides an explicit statement of the above in the parameter space {E ; h) and is proven in 
section ED 


Theorem 2.1. Let E > 1 and a > 0 be fixed. 

(i) Existence of positive solutions. There exists a unique h* = h*{E,a) such that 

- If h < h*, then there is no positive stationary solution and there is extinction of prey for the 
ODE system. 

- If h > h*{E,a), then there exist two positive solutions for the ODE system {u,v) and {u*,v*) 
with u < u*. 

(ii) Stability of the solutions. Let h > h*{E,a). The solution {u,v) is always unstable. 
Moreover, there exists a unique h**{E,a) > h*{E,a) such that 

- If h > h**{E,a) then {u*,v*) is stable. 

- If h*{E,a) < h < h**(E,a), the stability of {u*,v*) depends on r. It is unstable if r is small 
enough and stable otherwise. 


Remark 1. Our calculations show that the gap between h* and h** is very small, so that, roughly 
speaking, {u*,v*) is stable whenever it exists. However, if h belongs to the conditional stability zone, i.e. 
h G {h*,h**), and if r is very small, stability is lost and the system becomes excitable. This explains, in 
particular, the presence of pulses for small values of r when dealing with spatial interactions (see section 


0). 


The map {E,a) i— )• h*{E,a) has the following properties, as proved in section 5.1 


Properties 2.2. Let a > 0 be fixed. The map E i—)• h*(E,a) is increasing and one has the explicit 
limits : 

h*{E,a) = 2+ 2^1 +a. 

Let E >1 be fixed. The map a i—>■ h*{E, a) is increasing and one has the explicit limits : 


lim h*{E, a) = 
E-s-l 


1 + 0 if a < 1 
2 yW if a > 1 


lim 

E^-\-oo 


h*{E, 0) = ^(2E-1 + 2y/E{E-l)^ 


hi{E) and lim h*{E,a) = +oo. 

d^+oo 


Figure illustrates the maps h* and h** 


and the possible outcomes for system 
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Bistable Zone 


Conditional Bistable Zone 

(0,1) 

is unstable 

Monostable Zone 


1 encounter rate {E) 


Figure 1: Description of the dynamics of the ODE system @ in the E — h plane. If E < 1, the control solution 
(0,1) is unstable. In this zone there exists at least one positive steady state and the prey never 
disappears entirely. If E > 1, then the control solution (0,1) is always (locally) stable. Moreover, the 
E > 1 zone is the union of three subzones. Below the h* curve, (0,1) is the only non-negative steady 
state and is a global attractor : this is a monostable zone. Above the h* curve, there are two additional 
positive steady states, one of which is always unstable while the second, denoted {u*,v*), may be stable 
or unstable. Above the h** curve, (u*,v*) is always stable. In this subzone, the asymptotic behavior 
depends only on the initial conditions : this is a bistable zone. Between the h* and the h** curves, the 
stability of depends on other parameters : this is a conditional bistable zone. For illustrative 

purpose, the size of this last subzone has been considerably increased. 


2.2.2 Analysis of the PDE system 


We wish to identify the parameter conditions required to obtain prey extinction in the PDE system 
(§. A simple stability analysis shows that, if E < 1, then invasion occurs in the PDE system. If E > 1, 
then the situation for the PDE system is more complex. In this situation, the spatial structure and 
diffusion processes result in additional conditions for extinction. The rationale is explained in detail 
below. 

Let us assume that there is a positive stable stationary solution of 0 denoted by {u *, v* ) and that 
the initial condition u{x, 0) is close to u* at some places x and close to 0 at other places. Since both 
u* and 0 are stable, the demographic phenomena lead to an agregation near u* and an agregation near 
0. However, diffusion allows individuals to move around in space, so one of 0 or u* may be the final 
global attractor. In other words, there may be a (stable) traveling wave joining u* to 0. The direction 
of this wave, given by the sign of the speed of the wave, indicates whether extinction or invasion occurs. 
However, there are difficulties associated with this argument. 

- There can be no homogeneous stationary solution of 0, only stable heterogeneous positive sta¬ 
tionary solutions. In other words, it is possible that h < h*{E,a) without control occurring. 

- Even in the case of bistability {h > h*{E,a)), the bistable system 0 is neither competitive nor 
cooperative. Little theoretical knowledge is available concerning the occurrence of traveling waves 
in such systems, with even less known about the stability and direction of the wave. 

Using super and subsolutions, we show here how to obtain the conditions sufficient (but not neces¬ 
sary) for extinction and for invasion, based on well known scalar bistable PDEs. Roughly speaking, let 
(u(t, x),v{t, x)) be the solution of ([^. If we hnd a positive constant v such that, for any[^(t, x) G M"'' x M, 
v{t, x) > V then it comes 


dtu{t,x) — Axu{t,x) < u{t,x){l — u{t,x)) 


Eu{t, x) 

I -|- Ehu{t, x)~’ 


t > 0, X G M. 


2. It suffices that this condition occurs for t > to for some to > 0. 
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Let u be the solution of 


f dtu{t,x)-AMt,x) = u{t,x)il-u{t,x))-j^^^^^v, 

I u(0, x) > u{0, x). 

The comparison principle implies that u{t,x) > u{t,x). Now, if u{x,t) —)• 0 when t —)• +oo then 
u{t,x) —)• 0 when t — )• +oo and extinction of prey occurs (see figure [^(a)). Moreover, if u{x,t) = 
does not depend on x, then there is aspatial control (see figure [^(b)). 


Conversely, if we can identify a positive constant v such that v{x,t) < v, then it comes 

Eu(t x) 

dtu{t,x) — Axu{t,x) > u{t,x){l — u{t,x)) — -—;—^ > 0, x € 


1 + Ehu{t, x) 


Now, define u{t, x) as the solution of 


dtu{t, x) - A^u{t, x) = u{t, x)(l - u{t, x)) - 
u{0, x) < tt(0, x). 


( 6 ) 


The comparison principle implies that u{t, x) < u{t, x). It follows that if for some x G M, lim sup u{x, t) > 0, 

t —^ ~1“00 

then limsupu(t,x) > 0 and there is (non uniform) invasion (see figure!^, 
t—>-+oo U 

These arguments give rise to the following theorems yielding sufficient conditions, in terms of the 
parameters E, a and h, for extinction or invasion to occur. All theorems are proven in section We 
begin with a sufficient condition for uniform extinction. 


(a) Extinction 



(b) Uniform extinction 



Figure 2: (a) Sufficient condition for extinction. The solution u{t,x) is majored by a supersolution u{t,x). If 
lim u(t,x) = 0, then lim uit.x) = 0 and there is extinction. 

i-)-+oo t-)-+oo 

(b) Sufficient condition for uniform extinction. The solution u{t, x) is majored by a supersolution (/)(t) 
which does not depend on x. If lim ^(t) = 0, then lim u{t,x) = 0 uniformly in x and there is 

t—¥-\-oo t—>-+oo 

uniform extinction. 


Theorem 2.3. Let E > 1 and define 

hi{E) = 1 (2E - 1 + 2VE(E-1)) . 

If h < hi{E) then there is uniform extinction. In other words, for any initial condition verifying Q, 
there exists (j){t) > 0 such that any solution {u{t,x),v{t,x)) of (§ verifies 

Vx G M, 0 < u{t,x) < fi{t) and lim <f>{t) = 0. 

t^+oo 
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Invasion 



FIGURES: Sufficient condition for invasion. The solution u{t,x) is minored by a subsolution u(t,x). If 
linisup:u(t, x) > 0, then limsupw(t, a;) > 0 and there is invasion. 

>- + oo >- + oo 


If h > hi{E), then there can be invasion or extinction. The following theorem gives a sufficient 
condition for extinction to occur. 


Theorem 2.4. Let E > 1 be fixed and let {u,v) be the solution of Q-Q. Define v = l and let u be a 
solution of Q together with u{0,x) = u{0,x). There exists a unique h~ = h~{E) > hi{E), such that 

- If h < h~ (E), then Vx G M, lim u{t, x) = 0 

t^+oo 

- If h > h~{E), then Vx G M, lim x) = JI where Ji = Jl{E^ h) is a positive scalar. 

t^+oo 

As a consequence, if h < h~{E) there is extinction of prey. 


The map E h (E) verifies the following properties proved in section 


5.3 


Properties 2.5. The map E h (E) is increasing and admits the following explicit limits : 


lim h {E) = l ; liin h (E) = —. 

E^l E^+oo o 


Our last result gives a sufficient condition for invasion to occur. 

Theorem 2.6. Let E > 1 and a > 0 be fixed and let {u,v) be the solution of (§-(©. Define v = 
1 + o^ i^Eh u be a solution of Q together with u{0,x) = u{0,x). There exists a unique h'^ = 

h^{E,a) > h*{E,a) such that 

- If h < h'^{E, a), then Vx G M, lim u{t, x) = 0 

t^+oo 

- If h > h'^{E, a), then Vx G M, lim u{t, x) = /x where /x = fj.{E, h) is a positive scalar. 

X—>+oo — — — 

As a consequence, if h > h'^{E,a) there is invasion of prey. 

Finally, the following result specihes the behavior of the map /i"*". 

Properties 2.7. The maps E i—)• h~^{E, a) and a i—)• h'^{E, a) are increasing. For any E > 1, h'^{E, 0) = 
h~{E) and limQ_,.+oo h~^{E,a) = +oo. Finally, a > 0 being fixed, one has the explicit limit 


lim h~^(E,a) 

E^+oo 



Remark 2. The limit lim h~^{E,a) remains unknown. However, it can be proved that this limit exists 

E —^ 1 

and is greater than lim h*(E,a). 
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The results above are summarized in figure]^ In the domain {{E,h), E > 1, h~{E) < h < h~^{E,a)}, 
which we will refer to as the ‘transition zone’, it is not possible to draw any conclusions concerning 
whether prey invasion or extinction is likely to occur. Indeed, this zone can be separated into two 
subzones, according to the parameters values : 

Zone I = {{E, h), E > 1, max{h*{E, a), h~{E)) < h < h'^{E, a)}. 

Zone ll = {{E,h), E > 1, h-{E) <h< h*{E,a)}. 

In Zone I, our numerical simulations show non-monotonous traveling waves. In zone II, simulations show 
various types of behavior, including pulse and even heterogeneous positive stationary solutions. This 
phenomena are discussed in the section 



Figure 4: Description of the dynamic of the PDF system (H in the E — H plan. If i? < 1 there is always invasion. 

If i? > 1 there is a uniform extinction for h < hi{E), extinction for hi{E) < h < hr {E) and invasion 
for h'^(E,a) < h. The zone between h~ and h'^ is called the transition zone. This transition zone is 
splitted into two subzones : Zone I and Zone II, separated by the h* curve. In these two zones, both 
extinction or invasion of prey may occur due to various spatial phenomena. 
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3 Numerical study of the transition zone 


3.1 Influence of a. 

The mathematical results above demonstrate the influence of the parameters E and h, and, indirectly, 
that of the conversion rate a, on the long-term behavior of the system. More precisely, when E > 1, 
prey extinction or invasion may occur, depending on the value of h. Indeed, we can define two values 
h~ = h~{E) < = h~^{E,a) (see theorems 2.4 and 2.6). Extinction occurs if h < h~ and invasion 

occurs if h > h~^. When h G {h~, h~^) we observe richer dynamics, which may depend on other factors. 
We refer to this zone as the transition zone. Note that, as h~ is not dependent on a and h'^ is an 
increasing function of a (proposition 2.7), the size of this transition zone increases with increasing a. 


0 = 0 a — 0.5 




0=1 0=2 




Figure 5: Computation of h*{E, a), h~{E) and h^{E, a) for four values of a. The so-called transition zone lies 
between the red curve and the blue curve h~. The transition zone increases with increasing a. This 
transition zone can be split into two subzones separated by h* (black line). 


A first clue to the possible dynamics in the transition zone is provided by an understanding of the 
dynamics of the ODE system Q described in the theorem 2.1 The dynamic of Q is essentially de¬ 
pendent]^ on the position of h relative to h* = h*{E,a). When h < h*, there is no positive stationary 
solution, whereas for h > h* there are two positive stationary solutions, one of which, the larger of the 
two, is (nearly always) stable. 


The position of h* relative to h and h~^ provides a first description of the transition zone. By 


virtue of proposition 2.2, one gets the following. We always have h* < h~^ but the position of h* relative 
to h~ is dependent on a. On the one hand, from the facts that h*{l,a) > h~{l) for a > 0 and 
h*{E,0) = hi{E) < h~{E) for E > 1, we deduce that h~ > h* for large enough values of E and small 


3. The dynamics generally also depends on a quantity h**, defined in the theorem slightly greater than h* that is 
not taken into acconnt here for the sake of simplicity. 
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enough values of a. On the other hand, h* is an increasing function of a tending to +oo. We obtain 
that h*(E, a) > h~{E) for large values of a and any E > 1. 

Remark 3. When h~ > h*, which may occur for sufficiently small values of a, we see that taking space 
into account automatically increases the potential of extinction of preys. 

The transition zone can thus be separated into two subzones : one in which h < h* (Zone II) and one 
in which h > h* (Zone I). Figuresums up this discussion. As we will see below, both extinction and 
invasion are possible in each of these zones, but the phenomena at work differ considerably, according 
to whether h < h* or h > h*. These phenomena are studied in more detail below, using a numerical 
approach. 



encounter rate {E) 


Figure 6: Numerical computation of hcrit in the E — h space for the fixed values of a = 4, d = 1 and r = 1. 

The transition zone h S {h~, h'^) is split into two subzones separated by hern - Extinction occurs below 
hcrit and invasion occurs above. 


3.2 Extinction or invasion : influence of r and d. 

Our numerical analysis shows that both invasion and extinction are possible in the transition zone. 
When space is taken into account we see that h* does not separate the zone of invasion from that of 
extinction. These two zones are, indeed, separated by a new critical value of the handling time denoted 
hcrit £ {h~, h~^), which is dependent on E and a, of course, but also on the relative rates of growth (r) 
and diffusion (d) of the predator population ; 

hcrit — hcrit (.F ^ cx.fr ^ df 

As expected, when h > hcrit, prey invasion is observed, whereas extinction is observed when h < hcrit- 
Thus, higher values of hcrit ai'e associated with more effective predation and thus with less effective 
invasion by the prey. 

Figurecompletes the theoretical scheme represented in figure]^ by presenting an example of the 
curve {E, hcrit) for a particular selection of values for the parameters a, r and d in the E — h plan. 
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Like h* and h^, hcrit increases with both E and a. This naturally translates into the fact that, 
higher values of E increases the chance of meeting between predators and prey and that at higher a 
values, the predator is able to make greater use of the prey and can therefore eliminate it. 

Given the multiple dependence of hcrit on different parameters, hgure|^can only represent a parti¬ 
cular case, chosen for its simplicity. Figure shows the relationships between hcrit and d for fixed values 
of E and a and for various values of r. We see that hcrit increases with r. This translates into the fact 
that for small values of r, the predators growth rate is small and their effectiveness reduced. Conversly, 
hcrit (essentially) decreases with increasing d. This is due to the fact that for large d, predators spread 
into a zone in which prey are not present, resulting in a weakening predation. 
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diffusion rate (d) 

Figure 7: Numerical computation of hcrit with respect to d for four values of r with E = 2 and a = 4. For 
a given value of the parameters, there is invasion of prey if h > hcrit and extinction ii h < hcrit- 
Therefore, the larger hcrit is, the greater is the potential for extinction of prey. We see that hcrit is 
increasing in r and (essentially) decreasing in d. Thus, an increase of d or a decrease of r decrease the 
impact of the predation on invasive prey. 



Remark 4. For small values ofr and intermediate values ofd, predators may increase their effectiveness 
by increasing d. In that case, the predator growth rate being small, predators density remains high for a 
long time even if prey are absent. Now, if d is large enough but not large, predators may spread into a 
zone in which prey are not present and remain there at a high density and long enough to stop the prey 
to invade. This phenomenon may enable predators to form a barrier to prey’s movement, preventing 
thereby prey propagation. For too large values of d, the loss in predators effectiveness due to movement 
is too strong and the above phenomenon does not hold any more. This explains why, for small values of 
T, hcrit first increases and then decreases with increasing d. 
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Table 1: Summary of the four different situations involving space in the transition zone. Here, E = 2 and a = 4 
which yields h* « 5.4. The X axis represents space and the Y axis the concentration of species. The 
blue curve (top of each graph) represents the concentration of predators in space while the black curve 
is the concentration of prey. Extinction and Invasion Traveling Waves are abbreviated ETW and ITW 
respectively. 


3.3 Dynamics of the system in the transition zone 

We know that extinction occurs when h < hcrit, whereas invasion occurs when h > hcrit- It should be 
borne in mind that the existence or absence of non-trivial solutions that are homogeneous over space are 
dependent on the position of h with respect to h*. Consequently, the processes at work during extinction 
or invasion are highly dependent on these position. 

We will now describe the different dynamics occurring in the transition zone, summarized in table 
See section 2.2 for a precise definition of the various quantities described here. 


A] Invasion {h > hcrit) 

- Turing instabilities : h < h* and d ^ 1. As d increases, predators spread out, moving into areas 
from which the prey is absent, leading to a decrease in the size of the predator population. This 
phenomenon leads to a decrease in predator density throughout the space occupied by the predator, 
allowing the prey to survive in certain zones. We thus obtain a periodic distribution in space and 
a constant distribution over time of the densities of the prey and predator. Mathematically, this 
phenomenon is described by a Turing bifurcation. 

- Invasion traveling waves (ITW) : h > h*. The invasion is described simply by an invasion 
traveling wave : a wave of propagation linking the two stable solutions {u*,v*) and (0,1) in the 
direction of the positive solution (u*,v*). When r is large this traveling wave is monotone. By 
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contrast, when r <C 1 it displays a rich dynamics. When r <C 1 and d is not too large, we observe 
that hcrit is approximately equal to h* (see figure [^. This indicates that there is an invasion 
traveling wave if a positive solution exists. In this case, ahead of the front, v = 1 and, since r <C 1, 
the predator population increases very slowly. Besides, as h > h~ , the prey invades the space when 
the predator is at concentration 1. This leads to front advancing. Behind the front, the predator 
has had sufficient time to increase the size of its population and, therefore, to decrease the size of 
the prey population. As h > h* , the population of the prey decreases towards the positive solution 
u* and we observe a non-monotonous invasive traveling wave. 

B] Extinction {h < hcrit) 

- Pulse : h G [h~,h*] and r <C 1. Since r <C 1, the front of the wave is similar to the ITW 
described above. However, as h < h*, here is no homogeneous positive solution u*. The prey 
population therefore decreases to zero behind the front, whereas it continues to advance in ahead 
of the front. We thus obtain a pulse. 

- Extinction traveling wave (ETW) : h > h*. This corresponds to the simplest case described 
above. We observe a propagation wave linking the two stable solutions {u*,v*) and (0,1) in the 
direction of the control solution (0,1). 

Finally, figure presents the map in the E — h plane for fixed values of a, d and r. It furthermore 
specihes the possible dynamics in each zone. 



Figure 8: Different spatial dynamics in the transition zone. To ensure that all possible situations are represented, 
we choose a = 0.5, r = 0.01 and d = 10. Note that traveling waves of extinction (ETW) and traveling 
waves of invasion (ITW) may be obtained outside the transition zone. 


4 Conclusion and discussion 

4.1 Summary of the results 

We have shown that invasion occurs if E < 1. If predators do not encounter their prey they cannot 
control them. If E > 1, then extinction or invasion can occur, depending on the parameters h, E, a 
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and r. Uniform extinction occurs for h < hi{E) and extinction occurs for hi{E) < h < h~{E). Invasion 
occurs for h~^{E,a) < h. Thus, if h increases, we move from a zone of extinction without a considera¬ 
tion of space to a zone of extinction requiring a consideration of spatial aspects and then to a zone of 
invasion. For intermediate values of E, the zones of extinction increase with increasing E resulting in a 
higher potential of extinction, as hi,h*, h~ are increasing functions of E. When E is large, the zones of 
control do not depend on E any more because hi, h~, h~^ have hnite limits when E —)> +oo. Thus, E can 
only play a role in prey control if it takes intermediate values and if h is not too large. In summary, for 
low values oi E {E < 1) or high values of E or h, the outcome of the interaction (extinction or invasion) 
is independent of E. 

There is furthermore a transition zone splitted in two subzones, with various spatio-temporal phe¬ 
nomena and wherein both extinction and invasion can occur. The size of this transition zone greatly 
increases when the conversion rate a increases. Depending on the relative positions of these two zones 
with regard to the zones of extinction and of invasion, four spatial dynamics were identihed : extinction 
and invasion traveling waves, extinction pulse waves and heterogeneous stationary positive solutions of 
the Turing type. 


4.2 Biological interpretation of the main results 


We have shown that an increase in E increases the potential of extinction while an increase in h 
increases the potential of invasion. This translates the fact that a highly effective predator does have 
a high encounter rate and a small handling time. Furthermore, since h'^ is an increasing function of 
a, an increase in a decreases the potential of invasion and increases the size of the transition zone, 
which in turn increases the potential for the system to have complex dynamics. Finally, an increase of 
the diffusion rate d and a decrease of the amplitude of the predators growth rate r both increase the 
potential of invasion of prey. Thus, a generalist predator loses its effectiveness to exterminate invasive 
prey if it diffuses too fast or if it has a too slow dynamics. 


The above results are stated in term of adimensionalized parameters (see section 2.1). By choosing the 


appropriate spatio-temporal variables, we may dehne Du = ri = 1. Thus, the biological interpretations 
of d and r are accurate. Conversely, the dehnitions of the searching efficiency, E = EK 2 , the handling 
time h = and the conversion rate a = ^ complicate the biological interpretation of these three 
parameters. Thus, in addition to the above discussion about the influence of E and h, we now discuss 
our results in terms of the other biological variables : Ki, K 2 , 7 and r. The parameter h being increasing 
in the carrying capacity Ki of prey, prey with high carrying capacity show a high risk of being invasive. 
Conversly, h and E are respectively decreasing and increasing in the carrying capacity of predators i^ 2 - 
Thus, predators with high carrying capacity have a high potential to control prey invasion. Otherwise 
7 < 1 means that predator growth is mostly due to alternative prey while 7 > 1 implies that predator 
growth is due to consumption of the focal invasive prey. Finally, an increase in 7 and a decrease in the 
amplitude of the predator growth rate r yield an increase of a. Therefore, predators with a preference for 
the invasive prey or predators with a slow dynamics might display a complex dynamics. In particular, 
the likelihood of the system to exhibit a pulse wave is then important. 


4.3 The consideration of space often, but not always, increases the potential for 
control of pest invasion 

The model analyzed here was studied without taking space into account, except for a numerical 
exploration in the discussion, in an article by Magal et al. 2008 As explained in the introduction. 


models of identical structure have been proposed independently by Fagan et al. ||2002 and Chakraborty 


2015 We will now discuss our results in the context of these previous studies. Adding a spatial component 
to predator-prey systems makes any prediction about the controllability of the system difficult, as it 
then depends on the values of several parameters. The comparison between situations with and without 


15 









the consideration of space is epitomized by the distinction between h* , separating parameter regions 
of mono- and bistability in the ODE system, and hcrit, separating parameter regions of invasion and 
extinction in the PDE system. We will now focus on the case of E > 1, as values of E < 1 do not 
promote control, predators encountering prey too infrequently. 

If space is not taken into account, control occurs if /i < /i*, as 0 is a global attractor. This is still 
true in situations in which space is taken into account, if /i < hi where hi is smaller than h*. When h is 
between hi and h*, 0 is only a local attractor, so it is not possible to state that control is always attained. 
In this respect, adding consideration of space decreases the potential for control. Furthermore, when 
space is not taken into account, there is either extinction or invasion when h > h*{E,a), depending 
on initial conditions. Incorporating consideration of space changes the region where invasion occurs, 
for any values of the other parameters and for appropriate initial conditions, into h > h'^{E,a), with 
/i+ > h*. Thus, the consideration of space reduces the size of the zone wherein the invasion is certain 
and is detrimental to the invading prey. Finally, the relative levels of predator and prey diffusion also 
determine the potential for control. Our model shows that control is increased by predators being less 
mobile than prey. If predator mobility levels are too high, the predators become to thinly spread on the 
ground. For similar reasons, too high a level of prey mobility leaves the prey vulnerable to predators. This 
is entirely consistent with the experimental findings of Fagan et al. [2002 In conclusion, taking space 
into account can lead to an increase or a decrease in the controllability of invading prey by predators ; the 
addition of space to the model has no generic implication for considerations of predator-prey dynamics 
(see also Lam &: Ni ||2012 [Braverman et ah' 


2015). 


4.4 How can generalist predators reverse invasion by pest ? 

The originality of this study lies in its consideration of a generalist predator in a spatial context. 
When studying generalist predators, it is common practice to assume that the functional response is 


of type III, due to switching between prey species (Erbach et al. 2014; van Leuven et al. 2007, 2013 


Morozov & Petrovskii 2013). However, this approach is not mandatory, and other works (Basnet &: 
Mukhopadhyay [2014 ; Krivan fc Eisner [2006 Hoyle &: Bowers [2007 ) have considered a type H functional 
response. Altering our model to include a type HI functional response would be very costly in terms of 
understanding, because such responses lead to a loss of bistability. Its derivative would be null without 
prey, so some of our demonstration would fail and the analytical complexity would be greatly increased. 
However, traveling waves for specialist predator with type HI functional response are known to exist 0 
Sz Wu ||2008 ) which indicates that our result may be extended to this case. 

The complexity of analytical studies of spatial predator-prey interactions lies in the reaction terms 
being of alternative signs in the equations, making the study of the systems of equations essential 
1984b; [Huang et al. 


(Dunbar 


2003 Huang & Weng [2013). Other interactions, such as competition 


of two species (all negative) and symbiosis (all positive), are simpler, as their studies are similar to 

1994; [Alzahrani et al.' 


the study of a single equation (Volpert et al 


2012). This accounts for the 


slow scientific progress in this otherwise highly relevant topic. However, several major results have been 
obtained in recent decades, including those of the fundamental work of Owen Sz Lewis j2001 The finding 
of Owen and Lewis that predators can slow, stop, and even reverse invasion by their prey was based 
on the bistability of the prey-only dynamics of systems consisting of specialist predators attacking prey 
populations displaying Allee effects. By contrast, our work shows that the ability of generalist predators 
to control prey populations with logistic growth lies in the bistable dynamics of the coupled system. 
We also observe pseudo-Allee effects in our system, but their physics is quite different. An analysis of 
the ODE system identified parameter regions of monostable (extinction) and bistable (extinction or 
invasion) dynamics, but analysis of the associated PDE was able to distinguish different and additional 
regions of invasion and extinction. As a consequence, prey control was predicted to be possible when 
space was considered in additional situations other than those identihed without considering space. The 
reverse situation was also possible. None of these considerations apply to spatial predator-prey systems 
with specialist natural enemies. 
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5 Proofs 


5.1 Proof of theorem (2.1) 

Let E > 1 and a > 0 be fixed. For any /i > 0, the system (Q can be rewriten as 

= Qh{u)ifh{u) - v) 
mV = rv{gh{u) - v) 

where ©^(n) = /h(n) = ^(1 - n)(l + Ehu) and gh{u) = 1 + 


( 7 ) 


(i) Proof of the existence. Dehne H{h,u) = fh{u) — 9 h{u). For a given h > 0, a couple {u,v) is a positive 
stationary solution of (j^ if and only if u = fh{u) and 


u g] 0, 1[ is a solution of H{h,u) = 0. 


( 8 ) 


Now, fix u g] 0,1[. Since dhH{h,u) = m(1 — u) + a ^ i+Ehu ) ^ ^Tiat the map h i—)• H{h,u) is 

increasing. From FI > 1, we get H{0,u) < 0 and from u g] 0, 1[ we get lim E[{h,u) = +oo. The map 

h^+oo 

h —)• H{h,u) being continuous, this implies that for any u g] 0, 1[, there exists a unique h{u) > 0 such 
H{h,u) < 0 if /i < h{u), 
that < H{h,u) = 0 if h = h{u), 

H{h,u) > 0 if /i > h{u). 

The smooth function u i—>■ h{u) may be computed explicitly by noting that for any h > 0, the equation 
(l8|) is equivalent to the algebraic equation 


u g]0, 1[ is a solution of Ph{u) = 0 


wherein we have set 

Ph{u) = (1 — m)(1 + Ehu)^ — E{1 + Ehu + Eau). 
This yields the explicit formula 


h{u) 


Eu{l — u) 


^ ^ 

— + \/l + 4au{l - u)j + u - 1 


In particular 


lim h(u) = lim hiu) = + 00 . 

u^0+ u^l- 


This implies that the minimum of h{u) is obtained for some Ucrit £]0,1[. We define 


(9) 


( 10 ) 

( 11 ) 


h* = inf h{u) = min h{u) = h{ucrit)- (12) 

ms(o,i) ue(o,i) 

The definition of h* shows that if h < h*, then ([^ has zero solution. This also implies, together with 
the limits © and the continuity of tt 1 —)• h{u), that for any fixed h > h* the equation (© admits at 
least two solutions]^ (see the figure]^. In addition to this, for any hxed h > 0, one has deg(P/i) = 3 and 
Ph always admits a negative roots for P(0) = 1 — P < 0 and lim Ph{x) = + 00 . This implies that (1^ 

00 '—' 

admits at most two solutions]^ In conclusion, has exactly two positive solutions if h > h* and zero 
positive solution ifh<h*. This ends the proof of (i). 


4. Remark that h* > f;, because h < implies that /^ < 0 on (0,1) and ^ has no solution. 

5. These arguments also show that u i—>■ h{u) is decreasing on (0,u*) and increasing on (u*,l) ; for otherwise it is 
possible to choose h > 0 such that there is at least four different u £ (0,1) such that h = h{u), which is equivalent to Ph 
having at least 4 roots. See the figure 
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Figure 9: The four figures are computed ior E = a = 2, which gives h* « 4.36. The figure on the left represents 
the curve u i—>■ h{u) (in bold). Above the curve fh{u) > gniu) while below the curve /^(u) < gh{u). For 
a given h, the ordinate m of a point of this curve verifies fhiu) = gh{u) and corresponds to the positive 
stationary solution {u, fh{u)) of 0. The figures on the right represent the isoclines v = fh{u) and 
V = gh{u) for the three fixed values h = 3, h = 4.36 and h = 5. A positive stationary solution of the 
system Q corresponds to an intersection of these two isoclines. The system Q admits two positive 
solutions for h > h*, one (double) solution for the critical case h = h* and zero positive solution for 
h < h*. 


(a) Proof of the stability. Let h > h* he fixed. Let {u,v) be a positive stationary solution of 0. Since 
{u,v) verifies fh{u) = gh{u) = v, the Jacobian matrix at {u,v) reads 


J(n, v) 


®h{u)f'f^{u) -Qh{u) 
rvg'^{u) —rv 


hence 

det {J{u,v)) = rQh{u)v{g'f,{u) - f^iu)). 

From the proof of (i), we know that the system 0 admits exactly two positive solutions denoted 
respectively as {u, v) and {u *, v*) with u < Ucrit < u* and such that 


h = h{u) = h{u*), 


where u i—)• h{u) is given by (10). In particular, u and u* are the solution of 

H{h{u),u) = 0. 


Differentiating the equation (14) with respect to u gives 

duH{h{u),u) = —h'{u)dhH{h{u),u). 


(13) 

(14) 
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Thus, using the known fact that dhH{h{u),u) > 0, the identity (13) and the footnote]^ one gets 
duH{h,u) < 0 and duH{h,u*) > 0. Since u) =—( 5 ^(m)—/^(«)), this shows that det u)) <0 

and the instability of (u, v) follows. 

By contrast, one has det (J(u*, v*)) >0 and it appears that the stability of (u*,v*) is given by the sign 
of 

tr(j(u*,v*)) = Qh{u*)fh{u*) - rv*. (15) 

In order to highlight the dependence on h, for any h > h*, we note u* = u*{h) and we also define 
= tIt- From f^{u) = ([T^, we infer the following : 


~ If fj,{h) < u*{h) then {u*,v*) is asymptotically stable. 

- If fi{h) > u*{h) then the stability of {u*,v*) depends on r. More precisely, define Vcrit = ^ > 

0. (It is easy to show that Vcrit < !)■ 

- If r > Vcrit-, then {u*,v*) is asymptotically stable. 

- If r < Vcrit, then {u*,v*) is unstable. 

The sign of /u(/i) — u*{h) with respect to the parameter h remains to be found. 

On a first hand, the explicit expression of shows that /^* is decreasing and that = 0. 

Moreover, the definition (12) of h* yields u*{h*) = Ucrit and /^*(r(cHt) = g'h*{ucrit) > 0. Hence, g{h*) < 
u*{h). It is also clear that^ = lim fi{h) < lim u*{h) = 1. By continuity, we infer that the equation 

h —^~1“00 h —^~1“CXD 

li{h) = u*{h) has at least one solution in (/i*,+oo). 

On another hand, if //(/i) = u*{h) then H(h, fj,{h)) = 0, which may be rewritten as 

{Eh + 1)3 = dE{{Ehf + Eh{Ea + 1) - Ea). 

A simple analysis shows that this equation has exactly one negative solution, one solution in (O, and 
one solution in (^,+oo). Since h* > ^ (see footnote this implies that fi{h) = u*{h) has exactly one 
solution in {h*, +oo). We note this unique solution h** = h**{E, a). 

Finally, it is clear from the above arguments that /i(/i) < u*{h) if h € {h*,h**) and that g{h) > u*{h) 
if h > h**. This ends the proof of the Theorem. □ 


Proofs of properties 2.2 Similarly to the proof of the theorem 2.1, and to highlight the role of the 
parameters E and a, let us define 

H{E, a, h,u) = ^{1- u){l + Ehu) “ f 1 + « ^ ^Ehu 


From the proof of the theorem 2.1, we know that the quantity h* = h*{E,a) and the corresponding 


Ucrit = Ucrit{E, a) are characterized by the two equations 

H{E,a, h*, Ucrit) = 0 

duH{E,a,h*,Ucrit) = 0 . 


(16a) 

(16b) 


The implicit function theorem immediatly shows that the maps {E, a) i—?• {h*,Ucrit) belongs to ((1, Too) x 
[0,+cx));M^). 

- Proofs of the growth of Fi i— h*{E,a) and of a i—)• h*{E,a). 

Let a > 0 be fixed. Differentiate (16a) with respect to E and use (|16b|) yields 


OeH{E, a,h*, Ucrit) + dEh* {E, a) ■ dhH{E, a,h*, Ucrit) = 0. 
We already know that dhH < 0 and an explicit computation gives 


dEH{E^CX^h , Ucrit) — 


1 


> 0. 


E{1 + Eh*Ucrit) 

It follows that dEh*{E,a) > 0. Similar arguments show that dah*{E,a) > 0. 
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Computation of lim h*{E,a). 

E —>^1 

Let a > 0 be fixed. Since h*{-, a) is increasing and positive on (1, +oo), there exists a nonnegative 
scalar h*{a) such that h*{E, a) —)• h*{a) as E ^ 1. To compute this limit, denote P{E, a, h, u) = 
(1 + Ehu)H{E,a,h,u). From (16) and the definition of h*, one see that h* is the minimal value 
of h such that 

3m G [0,1], P{E^ a, /i, m) = duP{E, a, h, u) = 0. (17) 

In other words, h* is the minimal value of h such that P{E, a, h, ■) admits a multiple root in [0,1]. 
By passing to the limit —)• 1 in (|17[), we see that that h*{a) is the minimal value of h such that 


3m G [0,1], P(l, h, a, m) = duP{l, h, a, u) = 0. 

Explicit computations give 

P(l, a, h,u) = u + h{2 — h)u + a + 1 — h) . 

The multiplicity of the roots of P(l, a, h, •) need now to be discussed. 

• If /i < 2-y/a then 0 is the only root of P{1, a, h, •) and the multiplicity is one. 

• If /i > 2^/a then P{1, a, h, •) has two other real roots and may be explicitly written as 


(18) 


where 


P(l, a, h, m) = h ^ u{u — m_)(m — m+) 


u± = -^(ji-2± y^(2 - h)2 - 4(a + 1 - h 


- If h = 2y/a then m_ = m+ = 1 — and this multiple root belongs to [0,1) if and only if a > 1. 

- If 2^fa < h < a + 1 then the three roots 0, m_ and m+ are distinct and there is no multiple root. 

- If h = a + 1 then either 0 = m_ or 0 = m+, depending on the sign of a — 1. In both cases 0 is a 
multiple root. 

- Finally, if h > a + 1 then m_ < 0 < m+ and all the roots of P{f, a, h, •) have multiplicity one. 

The above discussion shows, using the characterization of h*{a), that 

7 \ f 1 + a if a < 1 

^ ’ \ 2y/a if a > 1. ^ ^ 

Note that a i—)■ h*{a) belongs to C'^([0,+oo),M"''). □ 

Computation of lim h*{E,a). 

E^+oo 

Since is positive and increasing, one has 

lim h*(E, a) = — 

E^+oo 

for some nonnegative number (wherein we have set g = +oo). 

However, since Ucrit{E, a) is bounded, up to a subsequence again denoted by E, one has Ucrit{E, a) — 
p for some positive number p (eventually depending on a). By passing to the limit in ( 16a| ), we 
deduce p > 0; for otherwise one obtains 0 = H{E,a,h* ,Ucrit) 3-oo. It follows that taking 
E —>■ +00 in (16b), one obtains /x = 1/2. Thus, by passing to the limit in (16a), one gets 

1 2 
£2 ~ £ ^ 
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and finally (£a being nonnegative) 


— — 2 + 2\/l + a; 


- Computation of h*{E,0). If a = 0 one gets 

P{E, h, 0, It) = —(1 + Ehu) iy—Ehu^ + u{Eh — 1) + 1 — S) . 


□ 

( 20 ) 


Standard computations show that P{E,h,0,-) admits two nonnegative roots if and only if /i > 
^{2E — 1 + y/E{E — 1) := hi{E) and one double nonnegative root if /i = hi{E). This shows that 


h*{E,0) = hi{E). □ 

Computation of lim h*(E,a). Since E i—)• h*{E,a) is increasing one gets for any E > 1, 

a^+oo _ 

h*{E,a) > h*{a). The explicit expression (fl^ of h*{a) implies lim h*{E,a) = +oo. □ 


5.2 Proof of the Theorem 12.31 

Let i? > 1 be fixed and {u{t,x),v{t,x)) be a solution of ([^ with initial condition (uo(x), uo(x)) 
verifying (©• 


Step 1. It is clear that u{t,x) > 0 for any t > 0 which implies 

dtv{t,x) — dAxv{t,x) > rv{t,x){l — v{t,x)), t > 0, x G M. 

Thus, using the comparison principle and the hypothesis that u(0,x) > 1, we deduce v{t,x) > 1 for any 
t > 0 and X G M. It follows that 


dtu{t, x) — Axu{t, x) < u{t, x)(l — u{t, x)) — 


Eu{t, x) 


1 + Ehu{t, X 

By the comparison principle, we infer that any solution u of 

Eu{t, x) 


-, t > 0, X G M. 


dtu{t, x) — Axu{t, x) = u{t, x)(l — u{t, x)) — 


I + Ehu{t, X 


t > 0, X G 


( 21 ) 


such that u(0, x) > uo{x) for any x G M satisfies 


Vt > 0, Vx G M, u{t,x) < u{t,x). 

In particular, let be the solution of the ordinary differential equation 


= ma - m) - 


t > 0 


together with the initial condition ())(0) = 1. 

^ is a homogeneous solution of (21) and from m(0,x) < 1 = (/>(0) we deduce 

Vt > 0, Vx G M, 0 < u{t,x) < 


( 22 ) 
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Step 2. Behavior of as t —)> +oo. 

It is clear that 0 is always a steady state of (22) and is asymptotically stable for E > 1. 

Now, let (/>o > 0 be a positive steady state of (22). <f>Q is a root of the polynomial P{E, h,0, ■) which 


is studied in the proof of the property |2.2| (see (20)). Hence, if h > hi{E) then (22) has two positive 
steady states that we denote as u~{E, h) < u'^{E,h) < 1 explicitly given by 


u^{E, h) 



(23) 


A linear analysis shows that for h > hi[E), u~{E, h) is unstable and W^iE, h) is (asymptotically) stable. 
Finally, classical arguments show that (/>(t) —)• u^{E,h) if /i > hi{E) and 4‘{t) —)• 0 if /i < hi{E). This 
ends the proof of the theorem. □ 


5.3 Proof of the Theorem 12.41 

Let FI > 1 be fixed and {u{t, x), v{t, x)) be a solution to ([^ with initial condition {uo{x), vo{x)) verifying 

(§. 

By the argument of the step 1 of the proof of |5.2| one already knows that : 

u{t,x) <u{t,x) (24) 


where u{t,x) is a solution of (21) with u{0,x) = uo{x). Moreover, from the step 2 of the proof of |5.2[ 
one knows that if h > hi{E) then the equation (HI) is bistable since ( |21| ) admits two stable nonnegative 
steady states : u = 0 and /x = u'^{E,h) < 1. We prove here that, in that case, there exists a trave¬ 
ling wave connecting Jx to 0 at a negative speed if and only if /i < h~{E) for some (implicit) number 
h~{E) > hi{E). This implies that for any x G M, u{t,x) —)■ 0 if h < h~{E) and u{t,x) —)• /x if /x > h~{E), 
which proves the theorem. 


Let h > hi{E) be fixed. It can be proven (see Fife [1979 , [Volpert et al. 1994 and the reference 
therein) that there exists a unique speed c = c{E, h) such that (21) admits a traveling solution of speed 
c which connects /x to 0. More precisely, there exists a profile U{^) verifying [/(—oo) = /x, t7(+oo) = 0 


and t/'(±oo) = 0, such that U{x — ct) = u{x,t) is a solution of (21). Moreover (see Fife 1979) this 
traveling wave describes the asymptotic behavior of all solutions provided the initial condition ([^ are 
verified. The theorem is proven by showing that there exists h~ (E) such that c{E, h) < 0 if h < h~ (E) 
and c{E, /x) > 0 if /i > h~{E). This result is a direct consequence of the two following lemma. The first 


lemma gives a characterization of the sign of c by an explicit function (see figure 10). 


Lemma 5.1. Define 


W{E,h,u)) = £ (s(l-s) 


Es 


1 -|- Ehs 


ds. 


Then sign(c{E, h)) = signiW{E, h, /x)) where /x := u^{E, h). 


Proof. Denote c = c{E,h) and let ^ = x — ct and f7(^) = u{x,t) with U{—oo) = /x, f7(-|-oo) = 0 and 
U'{±oo) = 0. We have 

-cU' = t/" + f/(l - f/) - = U" + H, U). 
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Multiplying by U' and integrating over M one gets 


=/r - T^) 

/■° dW 

= J -^{E,h,U)dU = -WiE,h,u+) 


and sign{c) = sign{W{E,h,u'^{E,h)) follows. 


□ 


□ 



Figure 10: Graph of u I— >■ —W{E,h,u). The steady states 0, u~ and correspond to critical points of the 
potential —W. A stable steady state is a local minimum of W and an unstable steady state is a local 
maximum. One sees that 0 and are two stable steady states. The sign of c characterizes which of 
them is the final global attractor : if c < 0, then 0 is the global attractor. If c > 0, then is the 
global attractor. 


The second lemma gives the sign of W with respect to h. 

Lemma 5.2. For any E > 1, there exists (a unique) h~{E) > hi{E) such that 

sign{W{E, /i, u'^{E, h))) = sign [h — h~{E)^ . 
Lemma 15.11 and 15.21 show that 


- There is invasion of prey for (21) {c{E,h) > 0) if /i G {h~{E),+oo). In that case, for any x G M, 
u{t, x) —>■ u'^{E, /i) as t —)• + 00 . 

- There is extinction of prey for (21) (c{E,h) < 0) if G [hi{E),h~{E)). In that case, for any 
X G M, u{t, x) —)• 0 as t —)■ + 00 . 

In particular, we infer from the inequality (24), that if G [hi{E),h~{E)), then for any x G M, 
u{t, x) —)• 0 as t —)> + 00 . This ends the proof of the theorem. □ 


It remains to prove the lemma 5.2 


5.1 


we know that 


Proof of lemma |5.2[ Dehne W{E,h) = W{E,h,u~^{E,h)). From the lemma 
sign{c) = sign{W{E,h)). We show here that there exists h~{E) > hi{E) such that sign{yV{E, h)) = 
sign{h — h~ {E)). 


Step 1. Differentiate W with respect to h gives 

dhW{E, h) = dhW{E, h, u+iE, h)) + duW{E, h, u+{E, h))dhU^{E, h). 
By the very definition of W and tt"*", one has duW{E, h,u'^{E, h)) = 0, which yields 

dhWiE, h) = dhW{E, h, u+{E, h)). 
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Explicit calculations show that 


2 1 
7 / li 7 / 

W{E, h,u) = ----- + —^ ln(l + Ehu). (25) 

Differentiate this expression with respect to h and denoting z = Ehu^{E, h) provides 

dhW{E,h,u^{E,h)) = - 21n(l + z) + • 

A standard analysis shows that the map z ^ z — 21n(l + z) + takes positive values for z > 0. It 
follows that the map h i— W{E, h) is increasing. 


Step 2. Recalling that u~^ is the largest roots of (20), one verifies that W^iE, /i) —>• 1 as /i —)• +oo and 
W{E, +oo) = 1/2 - 1/3 = 1/6 > 0. 


Step 3. This step consists in proving that for any E > l,W{E,hi{E)) = W{E,hi{E),u~^{E,hi{E))) := 
g{E) is negative. 


From the explicit expression (23) of u'^{E, h) and the definition of hi{E), we get 


u+{E,h,{E)) = ^{l- 


1 


= 1- E + y/E{E - 1). 


On a first hand, we have 

g'{E) = dEW{E,h{E),u+{E,h{E))) + dhW{E,hi{E),u+{E,h^{E))) • h[{E). 


Straightforward calculations give 


wherein we have set 


g'iE) 


(e/ii(E) 


2 


2 — ln(l + z) 



z = Ehi{E)u^{E, hi{E)) = (Ve + ^ {Ehi{E) - 1). 

A standard analysis shows that 2^^ < ln(l + z) for any z > 0. Thus g'{E) < 0 for any E > 1. 

On another hand, we have = 1 and 1) = 0, so that ^(l) = 0. It follows that g{E) < 0 for 

any E > 1. 


Conclusion For any F > 1, the map h i—)• W{E,h) is increasing and verifies W{E,hi{E)) < 0 and 
lim >V(F, h) = 1/6. By continuity, for any F > 1, there exists a unique h~ = h~{E) G (/ii(F), +oo) 

/l^ + OO _ 

such that sign(W{E, h,u~^)) = sign{h — h~{E)). This ends the proof of lemma 5.2 □ 


Proofs of properties |2.5 Let F > 1 be fixed. 
- Growth of 

Recall that h~{E) is characterized by 


lF(F,/i-(F),u+(F,/i-(F))) = 0 


(26) 


where the expressions of W and n"*" are respectively given in (25) and (23). By definition of W, 
one has 

duW (F,/i-(F),n+(F,/i-(F))) =0. 
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Hence 


dEW{E,h-{E),u+{E,h-{E))) + —{E)-dhW {E,h-{E),u+{E,h-{E)))=Q. 

The explicit computation of BeW and dhW are done in proofjT^and one gets 

BeW {E,h-{E),u^{E,h-{E))) < 0 and BhW {E,h-{E),u+{E,h-{E))) > 0 

so that 

^ BEW{E,h-{E),u+{E,h-{E))) 

dE^ ’ BhW{E,h-iE),u+{E,h-{E))) 

as needed. □ 

- Limit of h~{E) as i? —?■ 1. 

Let E > 1. One knows that h~{-) is increasing on (l,+oo) and minored by hi{E) > 0. Hence 
h~{E) admits a limit £~ as E ^ 1, and from hi{l) = 1, one obtains 

r > 1. 


The explicit expressions (23) of yield 


u {1,£ ) = 0andu^(l,f ) = 


l-£- 


Assume by contradiction that u'^{l,£ ) > 0. Then, for any h>l and 0 < u < u^{h, 1), we have 

1 


BuW{l, h, u) = u(l — u) — 


1 + hu 


> 0 and W (1, h, 0) = 0. 


It follows that ) = 0. Thus £ =1, which reads lim h (£’) = ! 

E^l 


This implies W(l,£ ,u~^(l,£ )) < 0, which contradicts (26). 

It follows that u~^(l,£~) = 0. 1 

Limit of h~{E) as E ^ +oo. 

One knowns that h~ [E] is inc 

(wherein we have set — = 0). By taking the limit E —>■ +oo in (23), we obtain 


□ 


One knowns that h (E) is increasing so that there is some £ > 0 such that lim h {E) = - 

E^+ca £ 


lim u^{E,h (E)) = ui := -^(1 + y/l — M) > 0. 

E^+oo 2 


and then, by taking the limit in (26), 




Thus, £ = ^ which reads, lim h (E) = 


E^+oo 


5.4 Invasion conditions, proof of Theorem 2.6 


As in the proof [5^ we start here by showing that 0 < u{t, x) < u{t, x) for some function u verifying 
a scalar reaction-diffusion equation (27) depending on E, h and a. Next, we show that there exists 
h^{E,a) such that limt_>.+oo m(L = fa > 0 when h > h^{E,a). This step is done using the proof of 
Theorem |2.4| together with an appropriate change of variables. 
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Step 1. Let > 1 be fixed. In general, using 0 < ri < 1, one has 1 <v <v with 

E 

V := v{E^ h,a) = 1 + a 


1 + Eh 


From the estimate v < v, we get u{t, x) < u{t, x) where u verifies 

Evu 


dtu = /S.xU + u{l — u) — 


1 + Ehu' 


t > 0, X G 


( 27 ) 


Denoting E = Ev and h = =, this equation reads simply 


dtu = AxU + u(l — u) — 


Eu 


1 + Ehu 


t > 0, X G 


(28) 


Note that, removing the , this equation is nothing but (21). As a consequence, the proof of the theorem 


2.4 implies the following result on (28). 


Lemma 5.3. Letu{t,x) he a solution of (28) verifying the initial condition ([^. There exists h~{E) > 
hi{E) such that 

- if h < h~{E) then for any x G M, lim u{t,x) = 0 (extinction), 

t^+oo 

- if h > h~{E) then for any x G M, lim u{t, x) = p(£^, h) > 0 (invasion). 

t^+oo 

Step 2. Recalling that E = Ev{E, h, a) and h = do depend on h, we set ^ = fJ.{E, h) = Ji{E, h) 

and the previous lemma gives the following implicit condition on h for invasion to occur. 

For any x G M, lim u{t, x) = //, provided : 

i —>-+00 — 

h > v{E,h,a)h~ {Ev{E,h,a)). (29) 

The following lemma gives an equivalent condition for this implicit condition to occur. This ends the 


proof of theorem 2.4 


Lemma 5.4. For any E > 1 and a>0, there exists (a unique) h'^{E,a) > h (E) such that : 

(29) holds true if and only if h > h~^{E,a). 


Proof of lemma |5.4| : Let E > 1 and a > 0 be fixed and define the function 

d^E,a{h) = h — v{E, h, a)h~ {Ev{E, h, a)). 


(30) 


By the construction of h via the implicit function theorem, one knows that EE,ai') is a function. 
Moreover, since h~ and v are bounded, we have EE,a{+oo) = +oo and 

J^E,a{0) = —v{E, 0, a)h~{Ev{E, 0, a)) = —(1 + aE)h~{E{1 + aE)) < 0. 

Therefore, it suffices to show that TE,a is increasing. 

One has 

J^E,aih) = 1 - dhv{E,h,a) (h~ (Ev) • 


From the expression of v, we infer dhV < 0 and from the proof of the properties 2.5 we know that 
> 0. It follows that Ef, ^{h) > 0. This ends the proof of the lemma. □ 
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Proofs of properties \2.7\ 

Proof of the growth of and h'^{-,a). 

h'^{E,a) is characterized by an equality in that is, 


J^E,a{h'^{E,a)) = 0 

where EE,a is defined in (30). Differentiating with respect to E gives 
dEh+{E, a) • E'E^^{h+iE, a)) + dE^E^h^iE, a)) = 0 

and then 

dEh+{E,a) ■ E'E^c.ih^{E,a)) = [{Ev + h-)dEV + v^] 
wherein we have set 


(31) 


dh- 

~dE 


(Ev). 


h =h {Ev), V = v{E,h~^{E,a),a) and Oev = dEv{E,h'^{E,a),a). 

One already knows that ^{h~^{E, a)) > 0 and that > 0. A direct computation shows that 

Oev > 0, which leads to dEh^{E,a) > 0 as needed. 


Similarly, differentiating (31) with respect to a gives, with obvious notations, 

dh- 


dah+{E, a) • E'E,c.ih+{E, a)) = d^v • + E 

and since daV > 0, one obtains dah'^{E, a) > 0. 

Limits of h^{E, a) as a —t- 0 and a —)• +oo. 

From v{E,h,0) = 1 we deduce h'^{E,0) = h~{E). Now, it is clear from the construction of 
that h'^{E,a) > h*{E,a). From the properties |2.2| we obtain h~^{E,a) —)> +00 as a —>■ + 00 . □ 

Limit of h~^{E,a) as FI —)• + 00 . 

First, recall that h~^{E,a) is characterized by (31), which reads 

E 


h'^iE, a) = ( 1 + a - ^ —r ) • h ( 

^ ^ V l + Eh+{E,a)J V 


h { E + a 


E^ 


1 + Eh+{E, a) 


(32) 


Since h^{-,a) is increasing, there exists £« > 0 such that liniE^+ooh'^{E,a) = Taking the 

limit E —)• +00 in (32) and using h~ (E) —>■ one obtains ^ ^ • (1 + aia)- The resolution of 

this equation ends the proof. □ 
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